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In this paper, a universal wave equation of the motion of microscopic particles in curved space-time is estab- 
lished, which satisfies the conservation of probability. Under the classical limit conditions, the wave equation 
returns to the dynamic equation of special relativity. The application of the wave equation in the evolution of 
a flat expanding universe has preliminarily proved the correctness of the equation, and the new wave equation 
may open a window for us to understand physics within the Planck scale. 


I. INTRODUCTION 


General relativity and quantum mechanics are two epoch-making contributions of modern physics, they have achieved brilliant 
success in their respective fields. The successful field of general relativity is the movement of matter in macroscopic large-scale 
space-time, in this field, quantum effects can be ignored. The successful field of quantum mechanics is the microcosmic world 
of atoms and subatoms, in this field, gravity can be completely ignored. However, in the events of strong gravity and micro-scale 
space-time, both general relativity and quantum mechanics will be involved. In order to describe such events, it is necessary to 
combine general relativity with quantum mechanics [1]. 

How to combine general relativity with quantum mechanics? One of these attempts is to quantize the gravitational field, such 
as the canonical quantum theory of gravity [2]. However, the quantized gravitational field cannot be renormalized, which is a 
divergent theory. Another attempt is to quantize space and time itself [3], which is a new view, but lacks experimental support. 

In this article, my goal is neither to quantize the gravitational field nor space-time, but to explore such a problem: in general 
relativity or curved space-time, what is the universal wave equation that the matter waves of microscopic particles satisfy? We 
require that the wave equation must meet the following conditions: 

(1) Ensure that the statistical interpretation of wave function is valid. That is, the square of the amplitude of the wave function 
is proportional to the probability of finding the particle. In other words, the wave describing the particle is a probability wave. 

(2) Guarantee the conservation of probability in curved space-time. 

(3) Ensure that the superposition principle of states is still valid. 

(4) The equation should contain the space-time metric tensor guv- 

(5)The equation should contain three basic physical constants at the same time: gravitational constant G, Planck constant h 
and speed of light c, indicating that it is a wave equation that combines general relativity and quantum mechanics. 

Obviously, neither Schrödinger equation nor Dirac equation can meet the above conditions at the same time. In this article, I 
will establish the wave equation of the motion of microscopic particles in the gravitational field(curved space-time), and lay the 
foundation for the further development of quantum mechanics. 


Il. ESTABLISHMENT OF WAVE EQUATION 


My idea is to extend Feynman’s path integral theory [4, 5] to curved space-time, and then derive the differential wave equation. 
Let y(x, yY, z, t) represent the probability amplitude of particle appearing at point (x, y, z) at time t, and W(x ; y, z, t) represent 
the probability amplitude of particle appearing at point (x ; y , z) at time ¢ , then there is the following relationship 


w(x, y, 2, t) = f Ke yesteal ye CU ye tda ay ae (1) 


where K(x, y, z, t; T, y j z : t) is called the propagator. According to Feynman’s path integral theory [4, 5], the propagator of 
particle from point A(z , y’, 2’) to point B(x, y, z) is 


K(B,A)=C X` exp{iS[r(t)]/h}, (2) 


all paths 
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where 


Sir(t)] = l TEDN (3) 


tA 


represents the action from point A to point B along the path r(t), L is the Lagrange of particle, C is the appropriate normalization 
constant. ñ in formula(2) is to make the phase factor dimensionless. The sum in formula(2) includes all possible paths from A 
to B. 

In curved space-time, the wave function ~)(x, y, z,t,7) is introduced to represent the probability amplitude of the particle 
appearing at space-time point (x, y, z, t) at time T. T is called proper time, (x, y, z, t) is the space-time coordinate, and there is 
a relationship between them dr? = —ds* = —g,,,dxdx”, g,,, is space-time metric tensor. By extending formula(1) to curved 
space-time, we get 


W(a,y,2,t,7) = f K(a,y,2,t,730 ,y,2,¢,7 We ,y,2,t,7 )V—gdt dx dy dz’, (4) 


where g = |g,» | is the determinant of the metric tensor, and , /—gd*x is the invariant volume element. K (a, y, 2,t,7; L,Y, z t, T) 
is the propagator in curved space-time, how to calculate it? I also extend Feynman path integral formula(2) to curved space-time. 
Let the propagator of particle from point A(x ,y ,z ,t ) to point B(x, y, z, t) be 


-o 5 anf > 


all paths 


where 


TB TR dxe dxf 
S[r(7)] = / dt = f Gap i at (6) 
TA 


represents the action from point A to point B along the path r(r), L = \/—gag de“ dz? is the Lagrange function of the particle, 


C is a appropriate normalization constant. lp = (Gh/c?)!/? = 1.61 x 10~3%cm, which is the Planck length. After introducing 
Planck length lp, the phase factor in formula(5) is dimensionless. Like formula(2), the sum in formula(5) includes all possible 
paths from A to B. 

Formulas (4) and (5) are the basic formulas that Feynman’s path integral theory is extended to curved space-time. From them, 
the differential wave equation of particle can be derived. 

We consider the state of particle (x, y, z, t, T +€) at time 7 +e, which has the following relationship with the state of particle 
W(x ,y,z,t,7) at time T 


+00 +00 +00 +00 fe Se et yor ro? 1 1 fogt 
veawetr+e) =f | J K(2,Y,Z,t,T +E; ,Y,Z,t,T)Y(£,Y,z,t,T)y/=gdt da dy dz. (7) 


co 


When £ — 0, in this infinitesimal proper time interval, according to formula(5), the propagator can be expressed as 
K(2,y,2,t,7 HEL, Y, zZ, tT) = C exp E |. (8) 


where C is a constant to be determined. Substituting formula(8) into formula(7), we obtain 


too too pto too ie E E oe T E ee 
wewatr+e) =f f / C exp Ea Way ,2,t,7T)/—g(x ,y ,z ,t )dt dz dy dz. (9) 
oo oo P 


Lett = t+ No, r =+ m, y = y+, z =z+ 73, then formula(9) can be written as 


+00 +00 +o +00 
p(z, Yz tT +E) = l- f. C exp E [vetm utm tm t+m:7) 


-f- Teese dni ns ie eden 


(10) 


When € tends to zero, it can be seen from the constraint condition dr? = — Judr” dx” that the main contribution of the 
integral comes from the region 79,1,2,3 ~ 0. Taking the integrand function in the above formula as the Taylor expansion of 


No. 71572, N3, we obtain 
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(11) 


In the € — 0, 70,771, 2,73 — 0 limit, we ignore all infinitesimal terms of higher order and equation(11) can be recombined 
N as 


+00 +00 +00 +00 1 4 82 
w(x, y, 2, t, T) +t D f. l- i n (x „Y, Zt, T) —g(£,Y, z, t) + 31 g2 [VgV] 


+518 5 [VEW] +; Lag [VT] + nb [v=o (12) 
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+Í: ` -|mn3 ale la -| N23 + i w(x, Y, 4, t,7) V —9(z, Y, at) dnodnıdnzdns. 


Formula(12) holds, then 


+00 +00 +oo +o 
/ f / CV —g(z, y, 2,t)dnodmdnzdn3 = 1 (13) 


is required. 
We have known that the integral formula 


+00 inè +00 in? +00 ine +00 ing 
a = 14 
L exp E -| f exp i, z nf exp H dn l- exp Pe | dna = (2nl pie)’, (14) 


so we can take the undetermined constant C as 


1 1 ing ini ing ing 
= A 15 
G vV 7g (2rl pie)? ied ES “or EE E 2l pe EE 2l pe ( ) 

We then use the integral formula 
+00 +00 +00 +00 
/ f / Cnydnodmdn2dns =0, u = 0,1,2,3. (16) 
+00 +oo +oo +oo 
Pe dd Crunvdnodm dnedng = 0, p #v a7 
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f f / Cridnodmdnzdns = wari (18) 


the equation(12) is converted into 


Ow ilpe [ 0? ə? 
penzt tete = yanat + eee EVT + aV) o 
oO oa 3? ) 
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which is the differential wave equation of the motion of microscopic particles in curved space-time(gravitational field). 
(1) This is a complex equation, and the statistical interpretation of the wave function is still valid in the curved space. For the 
normalized wave function, which satisfies 


fwen, z,t, T)? VZgd'z =1, (21) 


p = |w/? represents the probability density, and which conforms to the probability conservation equation(I will prove it in the 
next section). 

(2) Wave equation(20) is linear, that is, if ~%, and %2 are both solutions of the equation, then the linear superposition of 11 
and %2, ay, + biz, is also the solution of the equation, which conforms to the superposition principle of states in quantum 
mechanics. 

(3) In the natural unit system, lp = tpc = tp, the Planck time tp = 5.38 x 10-445, so another form of wave equation is 
to replace lp in equation(20) with tp. Equation(20), which contains three basic physical constants, gravitational constant G, 
Planck constant fi and speed of light c (because lp = (Gh/ evr 2), is wave equation of particle motion in general relativity. If 
Dirac equation combines special relativity with quantum mechanics, wave equation(20) is the combination of general relativity 
and quantum mechanics. 


Il. CONSERVATION OF PROBABILITY AND CLASSICAL LIMIT OF WAVE EQUATION 


We define the four-dimensional operator 


o o o 
= ér } z E 22 
ea w a my (22) 
where €;,, €y, €z, €; is the unit vector of four-dimensional rectangular coordinates, then 
3? 3? 3? O° 
be Ae oe oe 
so wave equation(20) can be rewritten as 
, Ow i2 2 1 
l — V : 24 
ipa NET (vV—g4) z” (24) 
By separating the amplitude and phase of the wave function, we let 
is 
Y= Retr, (25) 


where R and S are real numbers. Substituting formula(25) into equation(24), we get 


-àr 2S +20R-O9 + —=Dy-g: s); (26) 
OT 2 g 


əs ı 1 BR B 2 
— = —- (08) P 7 -OR. 
ap a) Fata R taa 9+ Fg 


(27) 


Equations(26) and (27) are completely equivalent to wave equation(24), now we will discuss their physical significance respec- 


tively. 


I first prove that equation(26) is the differential expression of the conservation of probability in curved space-time. The 


probability density p and flow density J are defined as 


p= |p|? = R, 
ilp : š 
J = E OY — y) 
Substituting formula(25) into equation(29), we get 
J = Ae, 


Both sides of equation(26) are multiplied by 2R, one can obtain 


Op 2J 
tLl-J4 =0 
OT Vg v59 
The definition of four-dimensional speed is 
da" 
Ue = 
dt 


From formula(30), the four-dimensional speed can also be expressed as 


U“ = (ge 


Equation(3 1) can be rewritten as 
Op ðJ” 2J" O 
+ + 
ðr Ox” ./—g Ox 


By integrating(34) with any volume element in curved space-time, we have 


(V—g) = 0 


[2v gdtx + | ZZ V-gd'x H [su (/—9)V—gdta = 0. 


The first term in formula(35) can be expressed as 
[Zz Op ae (fe zd | [eee 
V V = 


and the second and third terms in formula(35) can be combined as 
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(28) 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 
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which is the covariant divergence of the four-dimensional probability flow density J”. 
From formulas(36) and (37), equation(35) can be written as 
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The second and fourth items in the above formula are combined as 
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The reason for the last step in the derivation of the above formula is /=g = ./—g(2(rT), y(7), z(T),t(T)), and we have 


3 =g 2 3 =g dx $ O./—g dy + 3 =g dz P 3 =g dt 

OT Ox dr Oy dr Oz dr Ot dr 
Oy =9 at 
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So formula(39) is finally converted to 
o J pv =a] + | Juv =ada =0 
Using Gauss theorem in four-dimensional curved space-time, formula(42) can be written as 
2 J pv =a + $ V=gJ"d S, = 0, 


which is exactly the conservation equation of probability in curved space-time. 


(39) 


(40) 


(41) 


(42) 


(43) 


Now let’s discuss the classical limit of equation(27). Under the classical limit condition, Planck length lp — 0, the ZB terms 


in equation(27) can be ignored, then equation(27) can be simplified as 


ðS 1 


2 = 
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Because four-dimensional velocity U = OS, equation(44) can be rewritten as 
ðS ı 1 
ha SG, 
OT us 2 2 
Taking four-dimensional gradient on both sides of equation(45), and using U = US, we get 
OU 
=—+(U-O)U=0. 
OT +( ) 


Because the four-dimensional velocity can also be written as U = U(2(r), y(7), z(T), t(7), T), we have 


dU OU OUdx OUdy OUdz OUdt 


dr Or. be dr Oa dr Oe dr” Ot dr 
oU 
= -— +(U-DO)U 
OT 
At last, equation(46) can be expressed as 
dU 
— = 0 
dt , 


(44) 


(45) 


(46) 


(47) 


(48) 


which is exactly the equation of motion of free particles in special relativity. This means that when the planck length lp 
approaches zero, the gravitational and quantum effects disappear simultaneously, and the wave equation returns to the classical 


dynamics equation of special relativity. 


IV. THE EVOLUTION OF A FLAT EXPANDING UNIVERSE 
—A SIMPLE APPLICATION OF THE WAVE EQUATION 


Wave equation(20) is a complicated partial differential complex equation, which is difficult to solve in general, but can be 
solved in special cases. In this section, I will use the wave equation to discuss the dynamic evolution of the flat expanding 
universe. 

In a flat, expanding universe, the metric tensor is [6] 


-1 0 0 0 

0 a(t) 0 0 
m=! o o a(t) 0 oa 

0 0 0 a(t) 
which is called the Robertson-Walker(RW) metric, and a(t) is the scale factor of the expanding universe. 
Substituting RW metric into Einstein field equation, one can get 

34 = —4nG(p + 3p)a, (50) 
a = ae oa (51) 


p and p are the density and pressure of the medium, which are uniform in space and only change with time t, i. e. p = p(t), 
p = p(t). Equations(50) and (51) are called Friedmann equations, which are the basic equations of evolution of the universe. So 
we have two independent equations(50) and (51), but three unknown functions p(t), p(t) and a(t), thus they are incomplete. In 
the traditional discussions of the evolution of the universe, we need to artificially add a medium’s equation of state p = p(p) to 
form a complete set of dynamic equations. In a complete theory, the equation of state should be the natural result of the theory, so 
the theory which only contains Friedmann equations is incomplete. Now, we have an independent motion equation(20), because 
the wave equation is only applicable to the motion of non-zero mass matter, for the evolution stage of the universe dominated by 
non-zero mass matter, the Friedmann equations and the new equation of motion(20) together can form a complete equation set. 
In comoving coordinate system, dr = dt. To separate the amplitude and phase of the wave function, we define 
iS(t) 


w(t)=R(t)e  . (52) 


Substituting RW metric into equations of motion(26) and (27), we obtain 
. 1 p .. 6R . 
R= -3 RS+2RS + 7 (53) 
1 BB 3B., N 312. 3 
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Therefore we have a set of simultaneous equations composed of four independent equations(50), (51), (53) and (54). 


4nG 
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It is seemingly that we have five unknown functions p(t), p(t), a(t), R(t) and S(t), only four independent equations, which 
are still incomplete. However, considering the statistical interpretation of wave function p(t) = ||? = R? (t), there are only 
four independent unknown functions, so the set of equations (55) is complete and can be solved. 


The exact solution to equation system(55) is 


1 2 
a= (sex: + =) 
p(t) = 0 
a(t) = Ce (VérGt + C1) - (56) 
1 
RO) = Tae Cy 
S(t) =t+C; 


where C1, C2 and C3 are integral constants, which are determined by initial conditions. From the solution of the equations, 
we can see that the equation of state p = 0 and pa? =Const. are both natural corollaries of the complete dynamic evolution 
equation set(55). In addition, the four-dimensional velocity U? = Š = 1 and $? = U?Uo = —1 are also the natural results of 
the equations. 


V. EQUATIONS OF SPACE-TIME STRUCTURE AND PARTICLE MOTION 


In general relativity, particles move along the geodesic line in curved spacetime, and metric tensor g,» in the geodesic line is 
obtained by solving the Einstein field equation. So in general, we first assume that the distribution of matter is known (that is, we 
have known the energy momentum tensor T»), and then calculate the space-time metric tensor guy from Einstein field equation, 
and then calculate the motion of particles from the geodesic equation. Because the motion of particles is also a part of the energy 
momentum tensor, through the field equation the motion of particles will affect the structure of space-time g,,,, which in turn 
will affect the motion of particles again. In the most general case, the distribution of matter T,,,, and the space-time structure gy 
should be both unknown, in this case, how can we calculate the space-time metric tensor g,,, and the motion of particles at the 
same time? In general relativity, we can’t do it, the main reason is that the equation of motion, which is a corollary of the field 
equation, is not independent. That is, in general relativity, we have only one independent field equation, and it is impossible to 
solve g,,, and Ty simultaneously. But the motion of microscopic particles essentially follows the laws of quantum mechanics, 
now we have a new wave equation(20) of particle motion, which is a basic equation independent of Einstein field equation. 
The new equation of motion(20) combined with Einstein field equation, together with the statistical interpretation of the wave 
function, can form a set of equations to solve g,,, and T, simultaneously. 

Let the wave function of particle motion be 


p= Aet? , (57) 
where A and S satisfy the wave equation 
A 1 2A 
Of =) | arisen is 229-1008 |, (58) 
OT 2 J/g 
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then the energy momentum tensor of particle can be expressed as 


THY = purU” 
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When the wave function of a particle ~ changes, the energy momentum tensor T”” also changes according to equation(60). 
It is assumed that Einstein field equation is still valid during dynamic change, that is, if T#” changes at a certain time and space 
point, the metric tensor at the same time and space point will change immediately to ensure that the field equation is still valid. 
Substituting formula(60) into Einstein field equation, and combining motion equations(58) and (59), we obtain the simultaneous 


equations 
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where R*” is Ricci tensor, and the curvature scalar R = g” Ryy = Ri- The wave function of particle %(i. e., A and S) and 
the metric tensor g,,,, of space-time can be obtained simultaneously by solving the above simultaneous equations. The solution 
of the equations(61) is the real state of particle motion, that is, a state of equilibrium between the motion of particle and the 
space-time structure. 

(1). There are 12 unknown functions in the equation system(61): 10 metric tensors g,,,,, wave function’s amplitude A and its 
phase S. When the identity of curved space-time 


1 
(ne = gen) =0 (62) 
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is still valid, there are only 8 independent equations. The equations are underdetermined and can not completely determine 
Juv, A and S. According to the common method of general relativity, 4 conditions about g,,, can be added to form a complete 
set of equations together with 8 independent equations. These additional conditions are called gauge conditions or coordinate 
conditions. For example, a particularly convenient coordinate conditions are the harmonic coordinate conditions 
Tò = o"T* = 0. (63) 
(2). The new wave equation(20) contains a constant Planck length lp, when lp tends to zero, the wave equation returns to the 
classical equation of motion of special relativity, but within the Planck scale, lp cannot be ignored. Therefore, it is reasonable 
to believe that the new wave equation(20) is suitable for describing the motion of matter within the Planck scale. It should be a 
universal wave equation for particle motion in any space-time scale. If it is considered that the motion of matter always affects the 
corresponding changes of the space-time metric tensor through the Einstein field equation, the simultaneous equations(61) are 
still valid within the Planck scale. Within the Planck scale, the quantum fluctuations of space-time may make the identity(62) of 
curved space-time no longer hold, so the simultaneous equations(61) have 12 independent equations, 12 unknown functions, and 
the equations are complete. Therefore, the solution of the simultaneous equations(61) may open a window for us to understand 
the physics within the Planck scale, which needs further researches. 
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